Valley susceptibility of an interacting two-dimensional electron system 
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We report direct measurements of the valley susceptibility, the change of valley population in 
response to applied symmetry-breaking strain, in an AlAs two-dimensional electron system. As 
the two-dimensional density is reduced, the valley susceptibility dramatically increases relative to 
its band value, reflecting the system's strong electron-electron interaction. The increase has a 
remarkable resemblance to the enhancement of the spin susceptibility and establishes the analogy 
between the spin and valley degrees of freedom. 

PACS numbers: 71.70.-d, 71.70.Fk , 73.21.-b, 73.43.Qt 



Currently, there is considerable interest in controlled 
manipulation of electron spin in semiconductors. This 
interest partly stems from the technological potential of 
spintronics, namely, the use of carrier spin to realize novel 
electronic devices. More important, successful manipu- 
lation of spins could also impact the more exotic field 
of quantum computing since many of the current pro- 
posals envision spin as the quantum bit (qubit) of infor- 
mation [l], [2J, |3j . Here we describe measurements of an- 
other property of electrons, namely their valley degree of 
freedom, in a semiconductor where they occupy multiple 
conduction band minima (valleys) [Fig. [T^a)]. Specifi- 
cally, for a two-valley, two-dimensional electron system 
(2DES) in an AlAs quantum well, we have determined 
the "valley susceptibility", Xv, i-G., how the valley popu- 
lations respond to the application of symmetry-breaking 
strain. This is directly analogous to the spin suscepti- 
bility, Xsj which specifies how the spin populations re- 
spond to an applied magnetic field [Figs. 1(d)]. Our data 
show that Xv and Xs have strikingly similar behaviors, in- 
cluding an interaction-induced enhancement at low elec- 
tron densities. The results establish the general analogy 
between the spin and valley degrees of freedom, imply- 
ing the potential use of valleys in applications such as 
quantum computing. We also discuss the implications of 
our results for the controversial metal-insulator transition 
problem in 2D carrier systems. 

It is instructive to describe at the outset the expres- 
sions for the band values of spin and valley susceptibili- 
ties Xs,b and Xv,b 0- The spin susceptibility is defined 
as Xs,b — dAn/dB = gbpsp/Z, where An is the net 
spin imbalance, B is the applied magnetic field, gb is the 
band Lande g-factor, and p is the density of states at 
the Fermi level. Inserting the expression p = mb/nTi 2 for 
2D electrons, we have Xs.b = (jJ-B l^'Kh 1 )g\ > m\ } , where mb 
is the band effective mass. In analogy to spin, we can 
define valley susceptibility as Xv.b = dAn/de = pEij, 
= (l/nh 2 )mbE2,b, where An is the difference between 
the populations of the majority and minority valleys, e 
is strain, and i?2,& is the conduction band deformation 
potential ||. In a Fermi liquid picture, the interparticle 
interaction results in replacement of the parameters mj, 



<7b, and i?2,b |i| by their normalized values m* , g*, and 
E%. Note that Xs oc m*g* and Xv oc m*E^. 

Our experiments were performed on a high-mobility 
2DES confined to an 11 nm-thick, modulation-doped 
layer of AlAs grown by molecular beam epitaxy on a 
(001) GaAs substrate We studied two samples, each 
patterned in a standard Hall bar mesa aligned with the 
[100] crystal direction [Figs. 1(b) and (c)]. Using a metal 
gate deposited on the sample's surface we varied the 
2DES density n, between 2.5 and 9.5xlO n cm" 2 . The 
magneto-resistance measurements were performed in a 
liquid 3 He system with a base temperature of 0.3 K. 

In bulk AlAs electrons occupy three (six half) ellip- 
soidal conduction valleys at the six equivalent X-points 
of the Brillouin zone [Fig. OJa)]. For the 11 nm-wide 
quantum well used in our experiments, only the valleys 
with their major axes along [100] and [010] are occupied 
0; we refer to these as the Xand Y valleys, respectively. 
The application of symmetry-breaking strain along [100] 
and [010] splits the energies of the A" and Y valleys, trans- 
ferring charge from one valley to the other [Fig. Gib)]. To 
apply tunable strain we glued the sample to one side of a 
piezoelectric (piezo) stack actuator . The piezo polling 
direction is aligned along [100] as shown in Fig. GI C )- 
When bias Vp is applied to the piezo stack, it expands 
(shrinks) along [100] for Vp > (Vp < 0) and shrinks 
(expands) in the [010] direction. We have confirmed that 
this deformation is fully transmitted to the sample and, 
using metal strain gauges glued on the opposite side of 
piezo [Fig. QIc)], have measured its magnitude in both 
[100] and [010] directions HE|. 

We determined Xv via different techniques. The first 
involves measuring, at B = 0, sample's resistance (R xx ), 
along [100] as we induce tensile (compressive) strain 
and the electrons are transferred to the Y (X) valley 
[Fig- Gib)]. Because of the smaller (larger) effective mass 
of electrons in this valley along [100], R xx decreases (in- 
creases) as a function of Vp and saturates once all the 
electrons are transferred to the Y (X) valley We can 
use the onset of the saturation as a signature of the full 
electron transfer and, from the corresponding e, deter- 
mine Xv HH- Our second technique is based on measur- 
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FIG. 1: (a) The first Brillouin zone and conduction band constant energy surfaces for bulk AlAs. (b) Schematic diagram 
showing that electrons are transferred from the X to the Y valley upon the application of strain while the total 2DES density 
is independent of strain, (c) Experimental setup for the valley susceptibility measurements, (d) Energy diagram showing the 
spin (Ez) and valley {Ey) subband splittings with applied magnetic field or strain, respectively, (e) Sample resistance, R xx , 
vs the piezoelectric actuator bias, Vp at various filling factors v. The traces are offset for clarity. The dashed and dotted lines 
respectively mark Vp at which X and Y are equally populated, and Vp beyond which all the electrons are transferred to the Y 
valley. The top axis shows the strain as measured by the strain gauge, (f) Same data as in (e) but now plotted vs a normalized 
horizontal (top) axis ej ' B v , showing that the oscillations have the same period. The traces are offset for clarity and, for u—8, 10 
and 12, a smooth, parabolic background is subtracted from the traces to highlight the oscillations, (g) Plot of the coincidence 
index, i normalized by v, showing linearity with Vp. From the slopes of the lines fitted through the data points we obtain \v, 
and the i — intercepts give the value of Vp at which the X and Y valleys are equally populated (e = 0). Three data sets A, 
B, and C correspond to densities of: 2.8, 4.8 and 7.6xlO n cm -2 , respectively, (h) Schematic energy fan diagram showing the 
relevant energies: the cyclotron (Ec), Zeeman (Ez) and valley splitting (Ey)- 



ing the valley densities at a fixed Vp from the Shubnikov- 
de Haas oscillations of R xx as a function of perpendicular 
magnetic field, B±. The frequencies of these oscillations 
are proportional to the densities of the two valleys and 
provide a direct measure of Xv UJ3| • We quote data ob- 
tained from these two techniques later in the paper |l3j : 
here we would like to focus on a powerful variation of the 
second technique which also directly measures \ v . 



The technique is based on the "coincidence" of the 
2DES's energy levels at the Fermi energy (Ejr) and en- 
tails monitoring oscillations of R xx as a function of e 
at a fixed B± [14| . Examples of such data are shown in 
Fig.^e) for n =4.8xl0 11 cm~ 2 and various values of B± 
corresponding to even Landau level (LL) filling factors v 
— 4, 6,..., 12. Clear oscillations of R xx as a function 
of Vp are seen. These oscillations come about because 
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the strain-induced valley-splitting causes pairs of quan- 
tized energy levels of the 2DES in B± to cross at Ep. 
As schematically shown in the fan diagram of Fig. ^h) , 
there are three main energies in our system. The field B±_ 
quantizes the allowed energies into a set of LLs, separated 
by the cyclotron energy, Ec — heB±/m*. Because of the 
electron spin, each LL is further split into two levels, sep- 
arated by the Zeeman energy, Ez — [ipg*B. In the ab- 
sence of in-plane strain, each of these energy levels in our 
2DES should be two-fold degenerate. By straining the 
sample, we remove this degeneracy and introduce a third 
energy, the valley-splitting, Ey — e-Efi which increases 
linearly with strain, as illustrated in Fig. ^h) . When an 
integer number of the quantized energy levels of a 2DES 
are exactly filled, Ep falls in an energy gap separating 
adjacent levels (see, e.g. the shaded area of Fig. ^h) for 
the case of v = 8) and, at low temperatures, R xx exhibits 
a local minimum. As is evident in Fig.^h), however, at 
certain values of e, the energy levels corresponding to 
different valley- and spin-split LLs coincide at Ep. At 
such " coincidences" , the R xx minimum becomes weaker 
or disappears altogether |b|. From the fan diagram of 
Fig. ^h) we expect the weakening/strengthening of the 
R xx minimum to be a periodic function of e, or Vp, 
since it happens whenever Ey is an even multiple of Ec- 
Therefore we can directly measure Ey (in units of Ec) 
and determine the valley susceptibility: Xv — 4eB v /hAe, 
where Ae is the period of the R xx oscillations and B v is 
B± of the filling factor v at which the oscillations are 
measured. 

Before we discuss the measured values of func- 
tion of density, we point out several noteworthy features 
of Figs. 1(e)- (g) data. First, because of finite residual 
stress during the cooling of the sample and the piezo, we 
need a finite, cooldown-dependent Vp to attain the zero- 
strain condition in our experiments; this is about -250 V 
for the experiments of Figs. l(e)-(g) and is marked by a 
dashed vertical line in Fig. ^e). The data of Fig. ^e) 
themselves allow us to determine the zero-strain condi- 
tion: at Vp = -250 V, R xx at v = 4, 8, and 12 is at a 
(local) maximum while at v — 6 and 10 it is at a (local) 
minimum. This behavior is consistent with the energy 
diagram of Fig. ^h) when no strain is present. Second, 
the period of the R xx oscillations as a function of e is 
larger at smaller v (larger B±). In Fig.^f) we plot R xx 
as a function of e normalized by B v . It is clear in this 
plot that the oscillations are in-phase at v = 4, 8, and 
12, and that these are 180° out-of-phase with respect to 
the oscillations at v = 6 and 10. This is consistent with 
the simple fan diagram of Fig. ^h) (lif. To quantify 
this periodicity and also compare the data for different 
densities, in Fig. L^g) we show plots of the coincidence 
"index", i = Ey/Ec, divided by v, vs Vp; the indices 
are obtained from the positions of maxima and minima 
of the R xx oscillations and the index at the zero-strain 
condition is assigned to be zero. The plots in Fig. LUg) 
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FIG. 2: Valley and spin susceptibilities, normalized to their 
band values, vs density. \v was measured for two samples, SI 
and S2, and using various techniques, as indicated. \a data 
were measured in several AlAs 2DES samples, including S2, 
and are shown for the case when the X and Y valleys are 
equally occupied (open circles), and for when only one valley 
is occupied (open triangles). Lines are guides to the eye. 

show that the valley polarization is linear in e. This lin- 
earity means that Xv is independent of valley polarization 
or v (for v > 4) within our experimental uncertainty |l6| . 
Third, in Fig. L^e) R xx at all v shows a maximum at a 
particular value of e (dotted vertical line at Vp ~ 200 
V). This maximum marks the last "coincidence" of the 
energy levels. Beyond this value of e, the minority val- 
ley is completely depopulated and there are no further 
oscillations. We have confirmed this statement at lower 
sample densities where our available range of Vp allows 
us to go well beyond the valley depopulation e. 

Figure|2lsummarizes our main result, namely, the mea- 
sured Xv as a function of density. Xv is enhanced at 
all accessible densities and particularly at lower n. This 
enhancement is in sharp contrast to the non-interacting 
picture where Xv should maintain its band value at all 
electron densities, and betrays the strong influence of 
electron-electron interaction in this system. 

A comparison of the measured Xv with the spin suscep- 
tibility for similar samples |l7| supports this conjecture. 
As seen in Fig. [21 Xs and Xv show a remarkably similar 
enhancement. The increase of Xs as n is lowered is well 
established theoretically Il8| and exp erimentally for var- 
ious 2DESs [H H Lli HJ LH l23|. It occurs because 
of the increasing dominance of the Coulomb energy over 
the kinetic (Fermi) energy as n is decreased. To mini- 
mize its Coulomb energy, the 2DES tends to more easily 
align the electrons' spins so that they would spatially 
stay farther apart from each other. A similar mecha- 
nism is likely leading to the Xv enhancement. This is not 
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unexpected, if we regard valley and spin as similar de- 
grees of freedom. Indeed, in the standard effective mass 
approximation, where electrons near a semiconductor's 
conduction band minima are treated as quasi-particles 
with re-normalized (i.e., band values) of parameters such 
as (the effective) mass and g-factor but are otherwise 
free, spin and valley are equivalent, SU(2) symmetric de- 
grees of freedom. Interaction would then be expected 
to affect the spin and valley susceptibilities in a simi- 
lar fashion. Interaction-induced valley-splitting, akin to 
exchange-induced spin-splitting, is in fact, known to oc- 
cur in multi- valley 2DESs in the presence of B± [24| . Ev- 
idence was also recently reported for interaction-induced 
quantum Hall valley-Skyrmions, in direct analogy with 
spin-Skyrmions |25l |. 

Despite the similarities between the valley and spin 
degrees of freedom, there remains a puzzle. As shown in 
Fig. |2 within our experimental uncertainty, there is no 
spin polarization dependence of Xv, although there is a 
valley polarization dependence of Xs llZj- Our \v mea- 
surements are performed primarily at fixed v > 4 in a 
regime where the electrons are not fully spin polarized. 
But we have also measured Xv at v = 4 in tilted mag- 
netic fields, where the 2DES is fully spin polarized, and 
it agrees closely with Xv in the un-polarized regime (see 
"star" data point in Fig. EJ. However, as seen in Fig.0 
Xs in a single valley system is larger than in the two val- 
ley case |l7j ] . This disparity may point to a fundamental 
difference between the valley and spin degrees of free- 
dom, possibly reflecting the inadequacy of the effective 
mass approximation when one is dealing with interaction 
between electrons in different, anisotropic valleys located 
at different points of the Brillouin zone. 

We close by making some general remarks. First, uti- 
lizing the spin degree of freedom to make functional de- 
vices (spintronics), or as a qubit for quantum compu- 
tation, has been of much interest lately 0, 0, E3| The 
valley degree of freedom may provide an alternative for 
such applications. For example, it may be possible to 
encode quantum information into the "valley state" of 
a two-level quantum dot. This type of qubit could po- 
tentially have much longer coherence times than those 
of the spin-based GaAs qubits where the electron spin 
coherence times suffer from a coupling of the electron 
and nuclear spins (see, e.g., Ref. Q)- The results re- 
ported here demonstrate that the valley populations can 
be modified and monitored, and provide values for one 
of the most basic parameters of the system, namely its 
(valley) susceptibility. Second, the spin susceptibility of 
an interacting 2DES is expected to eventually diverge 
when n is lowered below a critical value at which the 
system attains a ferromagnetic ground state. There are 
also transitions to ferromagnetic quantum Hall states in 
samples with appropriate parameters. We believe that, 
under favorable conditions, e.g., at very small n or in the 
quantum Hall regime, "valley ferromagnetism" phenom- 



ena should also be prevalent. Finally, it is believed that 
there exists an intimate relationship between spin polar- 
ization and the 2D metal-insulator transition (MIT), and 
some have even suggested that the spin susceptibility di- 
verges at the zero-field MIT critical density [2(|. The 
similarity of spin and valley degrees of freedom suggests 
that the latter should play a role in the MIT problem 
as well. We have studied the MIT in the present sys- 
tem and our preliminary results indeed indicate that the 
2DES can be driven to an insulating phase by increasing 
its valley polarization. A full account of our findings will 
be reported elsewhere. 
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